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DIGITAL IMPLEMENTATION OF THE TF30-P-3 


TURBOFAN ENGINE CONTROL 
by David S. Cwynar and Peter G. Batterton 
Lewis Research Center 

SUMMARY 

The design of a digital control for the TF30-P-3 engine is described with de- 
tailed flow charts and programming methods . The control is designed to duplicate 
the existing hydromechanical control modes over the entire operating range . The 
control was implemented on a general purpose process control computer and required 
4833 words of storage and a 3.23-msec calculation time. Transient response of the 
digital control was evaluated by tests on a real time hybrid simulation of the TF30- 
P-3 engine. 

It is shown that the deadtime produced by the calculation time delay between 
sampling and final output is more significant to transient response than the effects 
associated with sampling rate alone . Transient response deterioration of the main 
fuel control limited update times to 75 msec for a 3. 23-msec calculation time, and to 
20 msec for a 20-msec calculation time. Extremely long (250 msec) calculation and 
update times could be tolerated without stability problems . 

Performance degradation of the exhaust nozzle control limited update time to 
50 msec for a 3.23-msec calculation time, and to 20 msec when calculation time and 
update time were identical. A delayed response of the afterburner light-off detector 
and exhaust nozzle overshoot with resulting fan oversuppression were the limiting 
factors . Extension of the update interval to 150 msec caused failure of the control 
due to a false indication by the blowout detector . 

Two methods of implementing the complex exhaust nozzle control are given , one 
of which is designed to reduce sampling interval requirements and allow a more 
lengthy calculation time . In addition a modification of the exhaust nozzle control to 
provide for either velocity or position servoactuation systems is discussed . 



INTRODUCTION 


Interest in digital control systems for airbreathing aircraft is increasing . New 
objectives of quieter engines, shorter takeoff and landing capabilities, higher effic- 
iencies with decreased mission costs and increased engine life cycles are placing 
higher demands on the control systems used. A close interaction of the airframe, 
inlet , and engine controls is becoming necessary , and the use of a digital computer 
in these flight systems is desirable as an efficient means of achieving this interac- 
tion. 

A starting point for a program utilizing the capabilities of digital control to ex- 
ercise improved control modes can be a basic digital program that reproduces the 
functions of the standard bill-of-materials (BOM) control modes for an engine . This 
computer program must permit efficient utilization of the digital computer's core ca- 
pacity and computing time so computing capacity is left for extended capabilities such 
as integrating inlet and engine controls , self-optimizing controls , and so forth . The 
basic BOM computer program must also provide control accuracy and dynamics com- 
parable to the hydromechanical system using only sampled information . 

The purpose of this report is to document a digital computer program which 
satisfies the previous objectives for the TF30-P-3 turbofan engine . This control du- 
plicates the function of the existing hydromechanical control, using control laws and 
logic paralleling that used by the hydromechanical control . The techniques used are 
sufficiently general that they may be used to produce digital versions of other , simi- 
lar controls with a minimum amount of effort . 

A description of the TF30-P-3 hydromechanical control is presented first . The 
digital problems of function generation and stability of internal dynamic loops are 
then considered as they apply to the TF30-P-3 control. A digital implementation of 
the control using a process control computer follows . This control is then used in 
conjuction with a real time simulation of the TF30-P-3 engine to evaluate the oper- 
ating characteristics of the finite state controller . Data obtained from the simulation 
is presented to illustrate the effects of calculation and update times and to compare 
the different exhaust nozzle programming schemes . 
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DESCRIPTION OF THE TF30-P-3 ENGINE AND CONTROL 


The TF30-P-3 is a twin-spool turbofan engine equipped with an afterburner . 
The engine, shown schematically in figure 1, includes a three-stage axial-flow fan 
mounted on the same shaft with a six-stage axial-flow low-pressure compressor . 
This unit is driven by a three-stage low-pressure turbine . A seven-stage axial- 
flow compressor driven by a single-stage air-cooled turbine makes up the high- 
pressure spool. 

The engine has a hydraulically actuated variable exhaust nozzle and seventh- 
stage (low-pressure) and twelfth-stage (high-pressure) compressor bleeds. The 
afterburner consists of a diffuser duct , a combustion chamber , a flameholder , and 
seven fuel spray rings arranged so as to identify five separate afterburning zones . 


Main Fuel Control 

A block diagram of the main fuel control is given in figure 2 . To help the read- 
er , the diagram is divided into three distinct areas separated by the dot-dash lines . 
All symbols used are defined in appendix A , The upper left portion of the diagram 
is the speed request calculation area . A desired high rotor speed Npt is developed 
as a function (MFC2) of power lever angle PLA, inlet total temperature Tt2 and 
pressure Pt2 (functions MFC 3 and MFC 4) . When the power lever reaches the 
afterburner request point of 70°, Npt becomes a function of Tt2 and Pt2 only. The 
desired speed Npt is upper and lower limited by the maximum and minimum speed 
schedules according to inlet total temperature and Mach number MN to form the re- 
quested speed NREQ . The maximum speed is increased somewhat by removal of the 
signal through KN, if afterburning is requested . 

The center portion of the diagram exhibits the acceleration and speed governor 
controls. Actual speed N2 is subtracted from NREQ , multiplied by a proportional 
gain Kg, and biased by the governor breakpoint fuel flow to burner pressure ratio 
Wf/Pb to form the proportional control’s desired Wf/Pb ratio POBL . This ratio is 
upper limited during acceleration and lower limited during deceleration to form the 
WFPB signal. The acceleration schedule MFC1 is a function of Tt2 and N2 and is 
designed to avoid turbine overtemperature and compressor stall, while the decelera- 
tion schedule is simply the constant Wf/Pb . A signal TL6 is generated for the ex- 
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haust nozzle control to indicate when WFPB is at least 9.87x10 (kg/hr)/ (N/m ) 
less than the acceleration schedule limit. 

The actual calculation of fuel flow is illustrated in the right-hand portion of 
figure 2 . The WFPB command is multiplied by burner pressure Pb to determine a 
fuel flow . Changes in this flow are rate limited , and the magnitude of the flow is both 
upper and lower limited . If no afterburner blowout signal (TL5 = 0) is received from 
the afterburner control and no engine shutdown is requested by a power lever posi- 
tion below 10° , this calculated flow becomes commanded fuel flow Wfc. During a de- 
tected afterburner blowout (TL5 = 1) , the desired Wf/Pb ratio is determined as a 
function of burner pressure by a blowout derichment schedule MFC 7 and becomes 
independent of speed. This modified ratio (WFBO) is multiplied by burner pressure 
and the control switches over to the resulting fuel flow command until TL5 is reset 
to zero . The TL5 signal may be reset by the exhaust nozzle area dropping below 
an area corresponding to 24° rotation of the exhaust nozzle pulley actuator, or a de- 
crease in speed below the afterburner turnon point (see fig. 3) . 


Exhaust Nozzle Control 

A block diagram of the exhaust nozzle control is presented in figure 3 . The di- 
agram is divided by the dot-dash lines into four distinct areas for easy reference . 

The basic function of the exhaust nozzle control is to maintain a desired burner 

pressure to turbine exit pressure ratio (Pb/Pt7) during afterburner operation . The 

s 

Pb/Pt7 schedule is a function of high rotor speed and burner pressure (ENC1 of 
fig . 2 and the Pb/Pt7 bias schedule) and will be called supression ratio Pb/Pt7 for 
reference . In addition the control performs afterburner light-off and blowout detec- 
tion , and provides a command to the afterburner fuel control . 

Logic to ensure proper sequencing and timing of the afterburner is also pro- 
vided by the exhaust nozzle control (center portion, fig. 3) . Prior to light-off, 
the blowout and light-off detectors, TL3 and TP2, the B/O & F/C relay valve lock, 
and the B/O & F/C relay valve are zero . When PLA is advanced beyond 70°, TL1 
energizes the P/L shuttle valve. When the engine reaches speed, a PR3 turn on sig- 
nal is received from the main fuel control . This signal will pass to the light-off de- 
tector and afterburner fuel control as T05 provided the engine is off the accelera- 
tion schedule by a small amount (TL6 = 1) . Once T05 has become energized , the 
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acceleration requirement is eliminated by the PXG relay valve operation. As soon 
as T05 becomes 1 , a small positive value (ERX) is passed to the zone fuel flow 
command piston (XOO) , This is illustrated in the right-hand portion of figure 3 . 

The ERX signal is integrated by the XOO piston to a small value determined by the 
Wf/Pb cam rise feedback schedule ENC6 and the commanded value XRQ from the 
correlation cam (ENC5) . This XOO position determines the initial light-off fuel 
flow. At the same time, IA14 goes high, arming the light-off detector (TP 2, center 
of diagram) to trip as soon as PPE increases 3 percent (L/O release schedule) above 
any minimum it may achieve before light off. The light-off detector will remain 
tripped as long as T05 remains high. 

When lightoff is detected (TP2 = 1) , the exhaust nozzle control is essentially a 
position servo with a proportional plus integral control on suppression ratio error 
PPE in the feedback path (lower left portion of fig . 3), The tripping of TP2 allows 
the integrator in the proportional plus integral control to integrate from its zero in- 
itial position and the loop around the exhaust nozzle actuator is closed , freeing the 
actuator from the hard closed command velocity which it had prior to lightoff . The 
output of the proportional plus integral control PIOT is a function of PPE, where 
PPE is calculated from the following equation: 

PPE = function ENC2 (B7E) 

where 



and K1 is a unity gain constant which decreases when B7E exceeds 15 percent. 

The commanded position into the feedback loop ARQ is either a function of power 
lever angle PLA , or an allowed position XAL derived from the zone fuel flow com- 
mand XOO , whichever is smaller . 

The exhaust nozzle position feedback AJP and the output from the proportional 
plus integral control PIOT are combined in the correlation cam ENC5. This cam es- 
sentially adds AJP and PIOT through variable , positive , or negative gains , and is 
used to produce the feedback to the actuator loop AFBX and to generate a command 
XRQ to the zone fuel control piston XOO . The primary function of this cam is to en- 
sure that the opening of the exhaust nozzle will lead increases in commanded after- 
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burner fuel flow . By generating a negative command to the XOO piston , this cam 
has the ability to cut back on afterburner fuel flow should PIOT become excessive . 

Whenever PPE drops below -15 percent, a blowout detection is possible. The 
blowout detector (left center portion of fig . 3) is armed by the exhaust nozzle open- 
ing 24° (TL4 = 1) , and the PR3 turnon command which powers T05. Therefore, a 
blowout may not be detected unless a lightoff was previously detected . All these con- 
ditions being met, TL3 will go high, flip the B/O & F/C relay valve to turn off T05 
and thereby reset the control for the next light-off detection . Note that if PLA is 
still advanced beyond 70°, the B/O & F/C relay valve lock will prevent the B/O & F/C 
relay valve from resetting, thus inhibiting future lightoffs until PLA drops below 
68°. The blowout signal is delayed 0,1 second and then passed to the main fuel con- 
trol as TL5 , 

In addition, a pop-open exhaust nozzle feature is provided for rapid thrust re- 
duction when landing (upper left portion of fig . 3) . It is activated by a squat switch , 
which detects when the plane is on the ground , and a power lever angle of less than 
22 °. 


Afterburner Fuel Control 

As seen from figure 4, afterburner fuel to burner pressure Wf/Pb ratios are 
scheduled for each zone as a function of the zone fuel flow command piston XOO . In 
addition, each zone’s fuel flow is set to zero unless XOO passes the zone cutoff point. 

Slightly beyond the cutoff point for zones 2 to 5 there is a hold position on the 
XOO piston (see lower right portion of fig . 3) . As soon as this hold position is 
reached , the TVL limit timer is activated to hold XOO at this position . A descrip- 
tion of the TVL limit timer operation is given in appendix B . This holding provides 
time for fuel to fill the zone fuel line , The amount of hold time is a function of burn- 
er pressure . If this hold timer were not implemented , the exhaust nozzle control 
would allow too rapid an increase in afterburner fuel flow since it would not "see” 
the large suppression ratio error caused by individual zone light-off transients until 
the deadtime needed to fill the fuel lines had elapsed . This additional phase lag 
in the loop would cause unstable operation; but after the line is filled , this phase lag 
no longer exists . 

The shape of the Wf/Pb piston feedback cam (ENC6 on fig. 3) is designed to 
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provide hysteresis in the XOO piston loop so that sharp decreases in XRQ due to 
light-off transients will not cause XOO to pull back . Also note that the Wf/Pb pis- 
ton velocity curve (ABC 8) is such so as to allow XOO to integrate rapidly in the 
forward direction , but much more slowly in the reverse direction . 


Bleed Door Control 


The seventh- stage bleed is controlled by inlet Mach number and is closed for 
Mach numbers below 1 . 7 (see fig . 5) . The twelfth-stage bleed operates on low com- 
pressor pressure ratio as follows: If the bleed is open, it will close when Ps3 > 


Kl^ieed ' + K2bieed an< * dPs3/dt >0.0, where Ps3 is the low compressor stat- 

ic discharge pressure and Kl^ lee( ^ and K 2 ^j ee( j are constants . If the bleed is 
closed, it will open when Ps3 < ^3^^^ • Pt2 + K4ki eec j and dPs3/dt <_ 0.0, where 
K3 bleed mA K4 bleed constants greater than Kl bUe(J and K2 bleed , respen- 
tively . The bleed will also open during a throttle chop (POBL < 12) or during an 
afterburner blowout (TL5 = 1) . 


CONSIDERATIONS IN DIGITAL CONTROL 
Function Generation 

One of the important aspects in building a control is the ability to generate the 
functions required . Functions of a single variable are easily generated by either 
linear interpolation between stored data points or by solving equations previously 
found by curve fitting the function . 

Functions of two or more variables are not so easily handled , however . One 
cannot assume that acceptable results may be obtained by simply keeping one varia- 
ble constant and generating two functions of the second variable to use in linear in- 
terpolation to obtain the desired result . Linear interpolation may be useless if the 
two single variable functions intersect or the first partial derivatives are discontinu- 
ous or contain points near infinity . 

An example of interpolation difficulty is depicted by the small portion of the ac- 
celeration schedule reproduced in figure 6. Here, discontinuous first partial deriv- 
atives are illustrated by the sharp breakpoints , and the single variable functions in- 


7 



tersect . By interpolating along lines parallel to the axes and between two functions 
of high rotor speed at constant Tt2, a third function representing a constant Tt2 
between the original two functions is obtained (dashed line) , Severe departure 
from the desired function occurs , partially because of the discontinuous partial der- 
ivatives , but mainly because linear interpolation forces all functions of constant Tt2 
between the original stored functions to intersect at the same point , For some func- 
tions , the errors introduced at the points of discontinuity or at a point where the 
partials are near infinity may be reduced by changing the interpolation axes or per- 
forming a linear transformation such as a change to corrected parameters . Refer- 
ence 6 offers one such solution. These techniques, however, will not correct errors 
introduced because of the intersection of the single variable functions . The solution 
to such problems can become quite cumbersome and time consuming , and may require 
the introduction of additional functions or variables . 

Fortunately the troublesome intervals which occur in the TF30-P-3 control 
either lie in a location where severe errors may be tolerated , or the errors which do 
result are so small as to be within acceptable limit's . In the case of the acceleration 
schedule , the troublesome areas occur at speeds below idle and , therefore , will only 
affect engine startup . In addition , the errors are such so as to produce a more con- 
servative acceleration schedule . The result is that the only detrimental effect of us- 
ing linear interpolation is a slight increase in startup time . The other schedules 
which produce similar difficulties can also be tolerated as producing either accept- 
able, or negligible discrepancies . These schedules are the maximum speed limit 
MFC 5 schedule, and the nonafterburning and afterburning biased speed request 
schedules of MFC 3 and MFC 4. 

Internal Dynamic Loops 

A second class of problems arising in digital control are those due to sampling, 
or to the discrete nature of the digital computer , Generally engine dynamics are 
slow enough such that control loops involving engine parameters are essentially un- 
affected by the digital computer's sampling limitations . Problems may arise if the 
control has self-contained fast response elements or control loops . Two such loops 
are identifiable within the TF30-P-3 hydromechanical afterburner control. 

Figure 7 depicts the result obtained when the afterburner control is linearized 
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and simplified using maximum gains for the nonlinear elements and logic states 
which produce worst case internal loop conditions . In this case the maximum values 
of functions ENC5, ENC6, ENC7, and ABC8 are defined as: 

max ENC5 (AJF) = Kac 
max ENC5 (PIOT) = Kpe 
max ENC6 (XQO) = Ke 
max ENG 7 (SVX) = Kav 
max ABC8 (EEX) = Kpv 

In creating this one-for-one digital model of the hydromechanical control , a sampler 
is required to "break up" any loop whose output depends upon either externally 
sampled or previously calculated values . Although samplers may be placed in many 
positions to accomplish the required result , analysis and calculation is simplified if 
a minimum number of samplers are chosen . For the case in point , placing a sampler 
after the feedback loop summing junction for loops A and B is sufficient . 

Zero-order holds are used on all analog outputs , and a deadtime of one sample in- 
terval is included wherever a value calculated during a previous update interval is 
used. 

Assuming the loops involving engine dynamics will remain stable under sam- 
pling , there are two identifiable loops which must be investigated for stable opera- 
tion. These loops are labeled as A and B in figure 7 . 

The stability of loops A and B can be investigated as follows . Referring to the 

AT. 

variables as labeled in figure 7 , the value of CR in the n u time interval was com- 
puted as 


CR (n) = [XRQ (n) - 1. 01 CR (n - 1)] Kpv • Kc • T + CR (n - 1) (1) 

where T is the update time in seconds. Using z-transform theory, equation (1) 
can be converted to the following z-transform: 

CR (z) = Kpv • Kc • T • XRQ (z) + (1 - 1. 01 Kpv • Kc « T) • [CR (z) . z _1 ] 


or 


CR (z) Kpv • Kc • T • z 

XRQ (z) ” z + 1. 01 Kpv • Kc • T - 1 


( 2 ) 
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For stability the root of the denominator of equation (2) must fall inside the unit cir- 
cle . The root is given by 


z = 1 - 1. 01 Kpv • Kc • T 


Hence, for stability 

-1 < 1 - 1.01 Kpv • Kc • T < 1 

Substituting the values of Kpv and Kc yields the result that loop A will be stable 
if 0 < T < 6 . 1 msec . 

For loop B analysis, the following relation (ref. 1) is used: 

AFBX (z) = G (z) 

ARQ (z) 1 + G (z) 

where G (z) is the z-transform of the network between SVX and AFBX . In evalu- 
ating G(z) a calculation deadtime must be added to account for the fact that the out- 
put does not occur immediately after sampling . Note that no deadtime was needed 
for the loop A calculation since no input or output was involved . Assuming a calcu- 
lation time equal to the update interval results in 

2. 2 Kav • Kb • Kac « T 
AJP(z) z (z - 1) 

ARQ (z) = 2. 2 Kav . Ks • Kac~ 

1 z (z - 1) 


Substituting the values of Kav , Ks , and Kac into equation (3) yields 

AJP (z) _ 1ST 
ARQ (z) z 2 - z + 18T 

The roots of the denominator of equation (4) are given by 



Restricting the roots to lie within the unit circle for T > 0 results in 

0 < T < 55. 5 msec 
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The previous results indicate that update times greater than 6 . 1 msec may pro- 
duce control instability . This limitation is severe and unnecessary . It should be 
pointed out that the problem has arisen because of the closed loop nature of the hy- 
dromechanical design . It is possible to avoid loops A and B by implementing a feed- 
forward equivalent of the control such as was done in reference 7 . However , it is 
not necessary to involve the complication of a redesign . One need only iterate around 
loop A a sufficient number of times to ensure stability . A sufficient number of times 
is simply some integer m such that T/m < 6 .1 msec . The procedure is to filter 
XRQ to produce m inputs to the loop for each update interval . For this control , a 
sample-and-hold filter was used . The m inputs thus produced were then used to 
calculate the loop as if it were operating at an update time of T/m seconds . Since 
the use of a sample-and-hold filter yields perfect knowledge of the current and fu- 
ture m - 1 inputs to the loop , the additional phase shift introduced by the filter can 
be effectively eliminated by advancing the iteration . For example , consider what 
happens if the current sampled value is held and put into the loop m times . The 
final result will be the same as if the update interval were T/m seconds, and the 
input assumed and held at the current value T(m - l)/m seconds earlier. This ad- 
vance will more than cancel the phase lag generated by the sample-and-hold filter if 
m > 2 . The additional phase lead helps minimize the phase lag produced by the 
sample-and-holds on the outputs . 

By using this iterative technique , one obtains an approximation to an advanced 
z-transform (ref 3.) of the loop involved; but unlike an advanced z-transform, the 
technique is applicable to nonlinear systems . The approximation becomes better as 
m is increased . This subinterval type of iteration is preferable to iteration tech- 
niques such as Newton-Raphson in that the calculation time is limited to a fixed num- 
ber of iterations . In Newton-Raphson , the error may not reduce to an acceptable 
level in the allotted number of iterations , causing unpredictable results . Through 
the use of subinterval iteration , then , internal dynamic loops need never restrict up- 
date times. 


DfGITAL IMPLEMENTATION OF THE TF30-P-3 CONTROL 


The techniques outlined previously were applied to implement the control from 
the block diagrams presented in figures 2 to 5 . The procedure is straightforward , 
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and a one-to-one correspondence exists between the hydromechanical and digital 
control modes . A description of the TVL limit timer implementation is given in 
appendix B , and a flow chart of the exhaust nozzle calculation (ENC subroutine) is 
given in appendix C . A detailed flow chart of the entire program is available from 
the author on request . 

To verify the digital implementation and to determine the effects of sampling , 
the control was programmed on a process control computer and used to run a real- 
time hybrid simulation of the engine , The configuration used is shown in figure 8 . 
Free stream pressure , temperature , and Mach number were entered into the simula- 
tion to set the operating condition . Mach number and update time were entered into 
the digital control from a teletype during program initialization. A complete descrip- 
tion of the hybrid engine simulation is given in reference 2 . 

Signals were received from the engine simulation as analog voltages . The con- 
trol outputs , in the form of analog voltages , were entered into an analog simulation 
of the control fuel lines and actuators . The actuator simulation is shown in figure 9 , 
and details of the digital computer are given in table I . The signal processing unit 
consists of buffer amplifers , readout devices , and patch boards used for convenience 
For the digital computer , approximately 4800 of the read/write memory cells were 
needed, and only one of the two multiplexer, sample-and-hold , digitizer units were 
used, limiting the maximum sample rate to 20 000 samples per second. Table II lists 
the various parts of the program , their approximate memory requirements and cal- 
culation time . 

An initialization program calculates the update time dependent constants used 
in the program . This initialization procedure is executed only once , and the con- 
stants calculated are stored as fixed values . An interval timer is used to issue an 
interrupt to the computer at a fixed rate equal to the specified update time . After 
receiving the interrupt , the computer begins sampling the input variables via a dir- 
ect memory access controller . The computer then idles for approximately 400 micro- 
seconds waiting for an interrupt from the direct memory access controller indicating 
that all variables have been digitized and stored . After receiving this second in- 
terrupt , control calculation is begun . As indicated in table II , approximately 
3.23 msec are needed to complete the calculation, leaving an additional idle period 
between the end of the control calculation and the next timer interrupt . During idle 
periods , the computer is available to display program variables to the operator . A 
timing diagram of the process is given in figure 10 . 
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OPERATING CHARACTERISTICS 


The control update interval determines how often the engine variables are 
sampled and how often the commands are updated . The control calculation time de- 
termines how much deadtime occurs between the sampling of the variables and the 
updating of the commands . Provided all internal dynamic loop problems have been 
properly handled , the maximum possible update intervals and calculation times will 
be determined by engine transient response . Since the engine is nonlinear , it would 
not suffice to evaluate these maximum s at a single operating point . For this reason, 
throttle bursts from idle to intermediate and from intermediate to maximum were used 
for evaluation purposes . All results were taken at sea-level static conditions . Unless 
indicated otherwise , a control calculation time of 3.23 msec was used . 


Main Fuel Control 

Figure 11 shows some important engine parameters during a throttle burst from 
idle to intermediate . Total mass flow is included as an indicator of engine thrust , 
since thrust itself was not available from the simulation . Inspection of the figure in- 
dicates that all parameters stay well within safe operating limits for update times as 
great as 249 msec . Update times beyond 249 msec could not be investigated , as this 
was the maximum time permitted by the scaling of the update time dependent param- 
eters . 

The important observation is that no significant performance change could be 
detected in the engine variables for update times to 75 msec . For the two rotor 
speeds , the response was nearly identical for update times to 150 msec . Such a broad 
range of operation indicates that digitizing the control has had negligible effect on its 
response to power lever commands. In addition, the major penalty for extremely long 
update times is an increased control response time . A verification of this fact can be 
seen from the high pressure compressor ratio response , shown in figure 12 . It is 
seen that stall margin improves as update time is increased . Although there is some 
ambiguity as to the exact location of the stall line , more than adequate stall margin 
exists for all operating conditions . 

Performance changes due to increasing update time show up best in the high 
pressure compressor ratio response . Observation of this response indicates a small 
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(2 percent or less) increase in stall margin for update times between 50 and 75 msec, 
this difference was too small to be conveniently included in figure 12 . With the in- 
creased sensitivity of this map , there is still no evidence of any benefit to updating 
faster than every 50 msec. 

Increasing the physical time it takes the computer to perform the control calcu- 
lation will introduce an additional deadtime into the control loop . This pure phase 
lag will have an effect on transient response in addition to those produced by in- 
creased update times . The response shown in figure 13 illustrates the point . Two 
responses to a throttle burst from idle to intermediate for a 150-msec update interval 
are compared to the 3.5-msec update, 3.23-msec calculation time response. The dif- 
ference between the two responses occurs in the amount of delay introduced by the 
calculation time . It is observed that the additional delay imposed by a 150-msec cal- 
culation time is greater than that delay produced by a 150-msec update interval 
alone. In general, this statement will be true since the phase lag of a sample-and- 
hold is equivalent to a deadtime of only half an update interval . The effect of making 
calculation time equal to update time is illustrated in figure 14 . For this case , per- 
formance degradation occurs at update intervals greater than 20 msec . 


Noniterative Exhaust Nozzle Control 

Since the linearized analysis of the internal dynamic loops of the exhaust nozzle 
control assumed worst case conditions , a test of a noniterative exhaust nozzle control 
was performed to determine where the practical limits would lie . Figures 15 , 16 , and 
17 illustrate the results. The control is seen to operate properly for a 3.5-msec up- 
date interval (fig. 15); but at a 10-msec interval (fig. 16), loop A begins to oscil- 
late. This is in reasonable agreement with the predicted value of 6.1 msec. The ef- 
fect of the oscillation is seen in the commanded afterburner fuel flow . At times this 
flow is proportional to XOO , and it is at these times that a loop A oscillation is re- 
flected in zone fuel flow . At an update time of 40 msec (fig . 17) , the oscillation be- 
comes so severe that it overides the PLA input command to the exhaust nozzle loop 
(loop B) , causing the nozzle to jitter. 


14 



Iterative Exhaust Nozzle Control 


The subinterval type of iteration as previously described was performed around 
loop A of the exhaust nozzle control to produce what will be referred to as the itera- 
tive exhaust nozzle control . It is this version which is represented by the exhaust 
nozzle control (ENC subroutine) flow chart in appendix C . An iteration time of 
5 msec was used . From the data to be presented , this time is seen to be sufficiently 
fast to obtain an acceptable approximation to the desired control , 

A throttle transient from intermediate to maximum is shown in figure 18 . At a 
50-msec update interval, performance of this control is still quite acceptable. The 
25-msec earlier start of zone 1 fuel flow is due to the calculation advance produced 
by the iteration around loop A. The advance is generated on the switching of T05 in 
a similar fashion as the advance which occurs on the XRQ input . This advance cal- 
culates to 45 msec, but only 25 msec is seen because of a potential lag of 50 msec on 
the sampling of throttle position demanding afterburner operation . The advance is 
desirable to offset this delay on PLA . The 25-msec early rise in exhaust nozzle 
position occurs because of an early light off due to the early fuel flow . 

The early light off of the remaining afterburner zones is partially due to phase 
lead generated by the iterative loop and partially due to the rounding which occurs 
when determining the required TVL limit timer hold time . The hold time must be an 
integer subdivision of the update time and is rounded down if the remainder is less 
than half of an update interval . 

Figure 19 gives a more complete set of engine parameters for a throttle burst 
from intermediate to maximum . At the 120-msec update time presented , there are 
several notable occurrences . First , there is no evidence of loop B becoming un- 
stable , even though the predicted maximum of 55 msec for a noniterative calculation 
has been exceeded . This discrepancy might be expected because of the fact that the 
actuator is held to its maximum slew limit during most of the transient , thus greatly 
reducing loop gain . A second consideration is that the calculation time is shorter 
than the assumed worst case value of one update interval , Zero calculation time in 
dicates instability wouldn’t occur until a 111-msec update interval was reached . The 
second feature of the response is that the exhaust nozzle starts opening late , even 
though zone 1 fuel flow begins early . This indicates the occurrence of a delayed 
light-off detection . Since the control is inhibited during the period between zone 1 
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turn on and light-off detection, the turn on of the remaining zones begins late. Like- 
wise , this delay has added to the phase delays produced by the longer update inter- 
val, aggravating an overshoot of exhaust nozzle position during the first 2 seconds 
of the transient . 

Light-off detection proves to be the ultimate restriction on update time . Fig- 
ure 20 illustrates what happens when update times of 150 msec are attempted. Be- 
cause of the phase shift added by the delay and a nearly unstable loop B , the con- 
trol becomes unstable in the time shortly following light-off. The sequence of events 
is as follows. Excessive overshoot of the nozzle causes severe oversupression . This 
oversupression causes an excessively negative supression ratio error such that the 
blowout detector is tripped . Once the blowout detector trips , the control turns off 
the afterburner and switches the main fuel calculation to the blowout derichment 
schedule . The rapid turnoff of the afterburner with the nozzle being open causes 
main burner pressure to drop . Fuel flow thus drops , with the process continuing 
until the exhaust nozzle closes sufficiently to reset the blowout detector . The normal 
fuel control then slowly brings the engine back to the prelight-off operating point . 

A practical limit on update time for the exhaust nozzle control is therefore de- 
fined as that point at which nozzle overshoot and oversupression become excessive . 
Practically , this point is around 50 msec . 

As was done for the main fuel control , the effects of increased calculation time 
on the afterburner control must also be considered . Figure 21 depicts these effects . 
The calculation advance produced by the iteration of loop A is cancelled by the de- 
lays of the increased calculation time . In addition , the assumed worst case calcula- 
tion time of one update interval used to determine the maximum allowable update time 
of 55 msec is now more accurate. Zone light offs are now increasingly delayed as the 
update and calculation times are increased . Acceptable operation can now only be 
obtained for a 20-msee update interval . If the iterative loop were to be replaced by 
a feed-forward type of implementation, phase lead such as that obtained by using an 
advanced z-transform , would have to be added in order to obtain proper operation 
at update times comparable to those of the main fuel control . 


Modified iterative Exhaust Nozzle Control 
For failsafe reasons , it may be desirable to operate the exhaust nozzle with an 
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external position type of actuation system . To demonstrate that the control can eas- 
ily accommodate such a system , the configuration of figure 22 was tested . In order 
to compensate for the inherent first-order lag produced by the position servo , a pro- 
portional plus integral control with AJV as its input was added . The net result to 
the control loop is a near return to the simple integration of AJV to produce posi- 
tion. To compare the two systems , refer to figure 7 . 

Since the actuator velocity limits were used to limit the AJV command, the 
integrator on AJV from the proportional plus integral control should be an accurate 
indication of exhaust nozzle position . Hence , two versions of the control were 
tried , one using actual position as the feedback , and one using the predicted posi- 
tion as the feedback . Note that the second system is a feed-forward control of the 
exhaust nozzle actuator . The performance of these two systems for a 50-msec update 
interval is seen in figures 23 and 24 . The lead on zone 1 turn on is due to causes sim- 
ilar to those of figure 18. The response of this control is nearly identical to that of the 
conventional iterative control and has acceptable performance at a 50-msec update in- 
terval , Figure 24 indicates that the predicted exhaust nozzle position is an accurate 
estimate of the actual position , and that the performance of the control using this pre- 
dicted value is also acceptable for a 50-msec update interval. 

The exhaust nozzle position servo actuation system has definite advantages . The 
accuracy of the predicted position provides a signal to compare with the measured 
position to determine if the actuator is functioning properly . This signal , along with 
the suppression ratio error , can then be used to effect a fail operational system . If 
both suppression ratio and the predicted versus actual position errors become large , 
the indication is that the exhaust nozzle is malfunctioning . A limitation against after- 
burner operation should then be imposed . However , if only the position error is ex- 
cessive , the indication is that only the position feedback transducer path is malfunct- 
ioning and feed-forward operation can then be employed for a restrictionless after- 
burner operation . The aforementioned fail operational system was not implemented 
but is offered as an indication of the versatility of the digital computer and the ad- 
vantages of employing a redundant electrical system . 


CONCLUDING REMARKS 

A digital control can easily duplicate standard hydromechanical turbojet control 
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functions. If fast response internal dynamic loops are required within the control, 
the problem of instability can be avoided by using an iteration technique which ef- 
fectively subdivides the update time into smaller increments. In addition, a method 
is available to determine the maximum update interval which ensures stability of 
these loops . 

A digital control similar in operation to the standard bill-of-materials control 
for the TF30-P-3 engine was proven feasible by having a process control computer 
implementation successfully control a real time simulation of the engine . The program 
required approximately 4800 words of storage and 3.23 msec for calculation. Using 
the programming techniques presented , calculation times up to 20 msec had negligi- 
ble effect on control operation and, for 3.23-msec calculation times, update times as 
large as 75 msec for the main fuel and 50 msec for the exhaust nozzle control could 
be tolerated before performance degradation took place . Sufficient computing capa- 
bility is therefore available for increased control complexity such as including inte- 
grated or self-optimizing, self-correcting control schemes. 

Lewis Research Center , 

National Aeronautics and Space Administration , 

Cleveland, Ohio, October 7, 1974, 

501-24, 
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APPENDIX A 


ABC1 - ABC5 

ABC6 

ABC 8 

ACR 

AFBX 

AJ 

AJP 

AJPC 

AJV 

ARM 

ARQ 

A$OT 

BL7 

BL12 

B7E 

CATCH 

CR 

d 

ENC 
ENC1 
ENC 2 
ENC 3 


SYMBOLS 

afterburning fuel flow ratio schedule for zones 1 to 5 , kg/hr 

Pb bias schedule 

Wf/Pb piston velocity schedule 

correlation cam output function of AJP and PIOT 

feedback from correlation cam to exhaust nozzle actuator loop 

2 

exhaust nozzle area, cm 

exhaust nozzle feedback pully position , deg 

exhaust nozzle feedback pully position command , deg 

velocity command to exhaust nozzle actuator, deg pully rotation/ 
sec 

light-off detector track on supression ratio error 
commanded exhaust nozzle pully position , deg 
power lever drive to exhaust nozzle 

low pressure compressor or seventh-stage bleed command 

high pressure compressor or twelfth-stage bleed 

supression ratio error 

TVL limit timer hold positions on XOO 

value of Wf/Pb feedback cam rise 

denominator 

subroutine which calculates exhaust nozzle control 
supression ratio schedule 
pressure ratio mixing cam characteristic 
proportional plus integral control integrator 
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ENC4 

ENC5 

ENC6 

ENC7 

ENC8 

ERX 

HOL 

HOLDTIME 

HOLT 

IA14 

IDLETRIM 

IPVD 

Kac 

Kav 

Kb 

Kc 

Kg 

KM 

KN 

KPB 

Kpc 

Kpv 

Ks 

K1 - K6 
LOOP 


schedule of TVL limit timer hold time 
correlation cam schedule 

afterburner Wf/Pb feedback cam rise schedule 

exhaust nozzle actuator velocity schedule 

supression ratio bias schedule 

error term in afterburner Wf/Pb piston loop 

index for TVL limit timer 

TVL limit timer hold time , sec 

number of update intervals in HOLDTIME 

logic for light-off relay 

addition to minimum speed schedule to trim for idle speed , rpm 

integral velocity in proportional plus integral control 

maximum correlation cam gain to exhaust nozzle actuator position 
signal 

maximum exhaust nozzle actuator velocity, deg/sec 
integral gain of bleed door actuators, percent/sec 
maximum gain of afterburner Wf/Pb feedback cam 
governor gain, (kg/hr)/ [(N/m )(rpm)] 
infinite gain 

maximum speed schedule bias , rpm 

bias on supression ratio schedule as function of Pb 

maximum correlation cam gain to PIOT 

maximum afterburner Wf/Pb piston velocity 

exhaust nozzle actuator integration rate, deg/ sec 

constants 

number of iterations per update interval on exhaust nozzle control 
implicit loop 
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LOPC 

MFC1 
MFC 2 
MFC 3 
MFC4 
MFC 5 
MFC 6 
MFC 7 

MILTRIM 

MN 

n 

NREQ 

Npt 

N1 

N2 

Pb 

PBT 

PICT 

PLA 

PLAx 

Po 

POBL 

PPE 

Ps3 


counter for exhaust nozzle control implicit loop (equals LOOP at 
start of each update interval) 

acceleration schedule 

speed request schedule 

nonafterburning bias schedule for MFC 2 

afterburning bias schedule for MFC 2 

maximum speed schedule 

minimum speed schedule 

blowout derichment schedule over range of Mach numbers (i.e . , 
function MFC 6 reduced to function of Tt2 only) 

trim on power lever for intermediate speed 

Mach number 

numerator 

speed request from power lever and schedules 
desired high rotor speed , rpm 
low rotor speed, rpm 
high rotor speed , rpm 

o 

burner pressure , N/m 

burner pressure biased to temperature 
output of proportional plus integral control 
power lever angle, deg 
power lever command, deg 
free stream static pressure 

2 

main fuel control desired Wf/Pb units, (kg/hr)/ (N/m ) 
modified supression ratio error 
static pressure , station 3 
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Pt2 


Pt7 

S 

SPR 

SVX 

T 

TCI - TC3 

TL1 

TL2 

TL3 

TL5 

TL6 

To 

TOl 

TO 2 

T03 

T05 

TP 2 

TSS 

Tt2 

TVLTIM 

T5 

Wf 

WFBO 

Wfc 

WFPB 


total pressure at engine face 

turbine exit total pressure 

Laplace variable 

power lever speed request 

exhaust nozzle actuator servovalve drive 

update interval, sec 

fuel valve time constants 

logic for power lever shuttle valve 

logic ofB/O & F/C relay valve lock 

logic of B/O arm valve 

blowout signal logic 

PXG turn on logic 

free stream total temperature 

PR 3 turn on logic 

logic of power lever shuttle valve 

logic of decel T/O arm valve 

logic of PXG relay valve 

logic of light-off detector 

error in exhaust nozzle control loop 

compressor inlet total temperature , °C 

subroutine which determines TVL limit timer operation 

turbine inlet temperature , °C 

main fuel flow , kg/ hr 

2 

blowout derichment schedule Wf/Pb units, (kg/ hr) / (N/m ) 
commanded main fuel flow , Kg/hr 

2 

commanded main fuel to burner pressure units , (kg/hr)/ (N/m ) 
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WFZ1 - WFZ5 


WSPA 

XAL 

XBL7 

XBL12 

XOO 

XOO 

XOOL 

XOS 

XRQ 

XWF 

XWFZ1 - XWFZ5 


afterburner fuel flow command for zones 1 to 5 , kg/hr 

2 

acceleration schedule Wf/Pb units, (kg/ hr)/ (N/m ) 

exhaust nozzle position allowed by zone fuel flow command posi- 
tion 

low pressure compressor bleed door position 
high pressure compressor bleed door position 
afterburner Wf/Pb piston position 

first time derivitive of XOO 

last value of XOO 

hold logic for TVL limit timer 

correlation cam drive to afterburner Wf/Pb piston 

main fuel flow , kg/hr 

afterburner fuel flow for zones 1 to 5 


Subscripts: 

bleed 

i 

m 

max 

min 

s 


pertaining to the bleed door control 

intermediate 

measured 

maximum 

minimum 

scheduled 
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APPENDIX B 


DESCRIPTION OF TVL LIMIT TIMER 
(TVLTIM SUBROUTINE) 

The TVL limit timer temporarily limits the travel of the afterburner Wf/Pb 
piston (XOO) at four catch positions. The purpose is to allow time for afterburner 
fuel flow to fill the afterburner fuel lines . The number of update intervals that XOO 
is held is determined by the negative number in HOL , where 

-HOLDTIME 


Prior to the TVL limit timer operation , the ENC subroutine calculates XOO 
using the last output XOS of the TVL limit timer . This calculated value of XOO 
along with the last value of XOO (XOOL) is used by the timer routine in the follow- 
ing manner: The timer is started when 

XOOL < CATCH and XOO > CATCH 

The timer is reset when the timer times out or when 

XOOL > CATCH and XOO < CATCH 

To start the timer , XOO is placed on the catch to its hold position , HOLT is set 
equal to HOL and XOS is set equal to -1. To reset the timer , HOLT is cleared to 
-1 and XOS is cleared to zero. During normal timer operation , the interval timer 
counter HOLT is incremented when leaving the implicit loop calculation of the ENC 
subroutine . If HOLT then becomes zero , the indication is that the timer has timed 
out during the last update interval and the timer is reset as indicated previously . 
The flow chart for Subroutine TVLTIM is given on the following pages . The catch 
positions are listed numerically next to their respective comparator logic elements . 
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TABLE I. - DIGITAL SYSTEM CAPABILITIES 


Digital computer 

Magnetic core memory size, words 

16 384 

Word length, bits plus parity 

16 

Memory cycle time, nsec 

750 

Add time, /usee 

1,5 

Subtract time, Msec 

1.5 

Multiply time, jusec 

4.5 

Divide time, jusec 

8.25 

Load time, jusec 

1.5 

Store time, jusec 

1.5 

Indirect addressing 

Infinite 

Indexing 

Total memory 

Priority interrupts 
Index registers : 

28 Separate levels 

Independent 

1 

In conjunction with lower accumulator 
Physical size, cm (in. ) 

1 

Width 

61 (24) 

Height 

157 (62) 

Depth 

76 (30) 

Interval timers 

Complement 

2 

Accuracy, clock pulses 

±1 

Clock rates, kHz 

572, 286, 160, 143, 80, 71.5, 
40, 35.75, 20, 10 

Counter 

16-Bit binary 

Output 

Priority interrupt to computer 

Analog acquisition unit 

Number of multiplexers, digitizers, 


sample- and- holds 

2 

Overall sample rate (maximum), kHz 

40 

Resolution of digital data, bits 

12 (plus sign) 

Output code 

Two's complement 

Number of channels 

64 

Input range, V full scale 

±10 

Input impedance, MS2 (shunted by 10 pF) 

10 

Maximum source resistance, Q 

1000 

Conversion time, jusec 

38 

Input setting time, jusec 

9 

Sample- and- hold aperture time, nsec 

500 

Safe input voltages, V 

±20 sustained 
±100 for less than 100 jusec 

Total error with calibration, percent 

0. 073 
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TABLE I. - Concluded. DIGITAL SYSTEM CAPABILITIES 


Frequency acquisition unit 

Number of channels 

10 

Nature of input 

Continuously varying or pulsatile 

Resolution of digital data, bits 

12 

Switch selectable clock rates, kHz 

20, 80, 100, 400, external 

Overall accuracy, bits 

±1 

Update rate 

Once per cycle of input frequency 

Maximum input frequency, kHz 

1 

Input amplitude range 

100 mV to 30 V (peak to peak) 

Analog output unit 

Total number of digital- to- analog 

26 

conversion channels 


Resolution (10 channels), bits 

12 (plus sign) 

Resolution (16 channels), bits 

11 (plus sign) 

Output voltage range, V full scale 

±10 

Output current (maximum), mA 

10 

Output impedance, Q 

<1 

Accuracy (12 bit), percent of full scale 

±0.1 

Accuracy (13 bit), percent of full scale 

±0.05 

Slew rate, V/ jU.se c 

1 

Settling time for 10- V step to within 0. 05 

20 

percent of final value, jusec 


Logical output unit 

Number of electronic switch outputs 

32 

Number of contact closure outputs 

32 

Maximum voltage, V 

30 

Maximum current, mA 

100 

Priority interrupt processor 

Number of channels 

10 

Input impedance, M2 

47 

Input voltage range, V 

±10 

Comparator switching 

Trigger on rise or fall 

Comparator hysteresis, mV 

Adjustable from 35 to 650 

Comparator output, V 

+7 

Monostable multivibrator: 


Pulse width, Msec 

0.3 

Pulse height, V 

+7 
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TABLE H. - CONTROL MEMORY AND TIMING REQUIREMENTS 


[Control calculation plus sample time, 3. 23 msec. ] 


Routine 

name 

Routine function 

Memory 

words 

reguired 

Memory cycles required 

Computation time, msec a 

Minimum 

Maximum 

Minimum 

Maximum 

FG 

Single variable function 
generation 

26 

64 

64 

0,048 

0. 048 

FG2D 

Two variable 
function generation 

79 

133 + 2 FG 

b 149 + 2 FG 

0.100 + 2 FG 

0. 112 + 2 FG 

TIME 

Timer and data 
acquisition control 

13 

22 

22 

0.017 

0. 017 

EOB 

Input scaling, calling 
of subroutines, and 
control outputs 

272 

236 

246 

0. 177 

0. 185 

SPREQ 

Speed biasing and limiting 
for inlet conditions and 
afterburner operation 

121 

67 + FG+2 FG2D 

69 + FG + 2 FG2D 

0.050 

0.052 

WFCOMP 

Main fuel calculation 
and supression ratio 
scheduling 

248 

272 + 2 FG + 2 FG2D 

316 + 4 FG + 2 FG2D 

0. 204 

0.237 

IPCTRL ! 

Exhaust nozzle integral pro- 
portional control, light- off 
detector, afterburner zone 1 
turn on 

84 

86 

113 + FG 

0. 065 

0.085 

ENC 

Correlation between after- 
burner fuel and exhaust 
nozzle 

180 

94 + FG + FG2D + TVLTIM 

b 121 + 2 FG + FG2D + TVLTIM 

0.071 

0.091 

TVLTIM 

TVL limit timer 

68 

21 

40 

0.016 

0.030 

BLEED 

Bleed door control 

56 

51 

61 

0.038 

0.046 

AB 

Afterburner fuel flow cal- 
culation 

272 

229 

b 825 

0.172 

0.619 

MAIN 

Program initialization 

238 

(c) 

(c) 

(c) 

<c) 

BLOCK 

DATA 

Data for schedules 

3176 

(d) 

(d) 

(d) 

(d) 

Total control calculation 

4833 

2639 

3774 

1.98 

2.83 


a Based on memory cycle time ol 750 nsec. 

b Based on independent variable of function changing only one tabulated value per update interval (see text). 
c Executed only Once upon program start. 

^Nonexecutable data. 
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High pressure 



Figure 1. - Station identification. 




Figure 2. - Main fuel control 













































































































Figure 4. - Afterburner fuel control block diagram for TF30-P-3 engine. 
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Figure 5. - Bleed door control block diagram. 
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Figure 6. - Inaccuracy produced by using linear interpolation to generate 
functions of two variables. (Not to scale. ) 



Figure 7. - Linearized model of exhaust nozzle control showing internal dynamic loops. 
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Update interval 


Calculation time- 


Timer Sampling com- 

mterrupt plete interrupt 


Reinitialize 
timer and 
start measure- 
ment sampling 


Initialize 
for implicit 
loop and TVL 
limit timer 
calculation 


Idle 

during - 
sampling 


Perform control 
calculation 


Output 

control 

parameters 


Idle until 
next timer — 
interrupt 


Time, 0 
Msec 


6 17 417 -3176 varies- 3230 t 

maximum time 
indicated 


Figure 10. - Timing diagram for typical update interval. When calculation time equals update interval, previous 
interval's calculation is output in 6- to 17 -psec interval, extending this interval by 54 Msec and the last two 
subintervals are eliminated. ) 
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Figure 11. - Effects of update interval on simulated engine response 



High compressor corrected weight flow 

Figure 12. - Effects of update interval on high compressor pressure ratio for throttle steps and chops between idle and 
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Figure 14. - Effects of extended calculation time and update interval on simulated engine response for th rottle burst from idle to military at sea-level static conditions, 
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Figure 16. -Throttle burst from military to maximum for noniterative exhaust nozzle control using 10-msec 
update interval. 


49 





Time, sec 

Figure 18. - Throttle burst from military to maximum for iterative exhaust nozzle control with acceptable control at 50-msec update interval. 
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Conventional control operating at 3. 5-msec update interval 
Modified iterative control operating at 50-msec update interval 
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Figure 22. - Block diagram of modified portion of exhaust nozzle control. 
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Figure 21. - Effects of increased calculation and update times on iterative exhaust nozzle control for throttle bursts from military to maximum. 
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